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Abstract 

Taylor's law (TL), which states that variance in population density is related to 
mean density via a power law, and density-mass allometry, which states that 
mean density is related to body mass via a power law, are two of the most 
widely observed patterns in ecology. Combining these two laws predicts that 
the variance in density is related to body mass via a power law (variance-mass 
allometry). Marine size spectra are known to exhibit density-mass allometry, 
but variance-mass allometry has not been investigated. We show that variance 
and body mass in unexploited size spectrum models are related by a power law, 
and that this leads to TL with an exponent slightly <2. These simulated rela- 
tionships are disrupted less by balanced harvesting, in which fishing effort is 
spread across a wide range of body sizes, than by size-at-entry fishing, in which 
only fish above a certain size may legally be caught. 
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Introduction 

Two widely observed empirical patterns in ecology are Tay- 
lor's law (TL; Taylor 1961; Eisler et al. 2008) and density- 
mass allometry (Blackburn and Gaston 1999; Jennings et al. 
2007; Belgrano and Reiss 2011). TL asserts that, in an ensem- 
ble of populations, the variance of the population density is 
a power-law function of the mean density of those popula- 
tions. Density-mass allometry asserts that population density 
is a power-law function of mean body mass, and can refer to 
single or mixed species as well as to individuals regardless of 
species. Composing TL with density-mass allometry predicts 
that the variance of population density should be a power- 
law function of mean body size and that the parameters of 
that power law should be predictable from the parameters of 
TL and of density-mass allometry (Marquet et al. 2005 and, 
independently, Cohen et al. 2012). Cohen et al. (2012) con- 
firmed this relationship using detailed forestry data. 

Taylor's law has been confirmed for hundreds of spe- 
cies or groups of related species in field observations and 



laboratory experiments (Reed and Hobbs 2004; Benton 
and Beckerman 2005; Marquet et al. 2005; Eisler et al. 
2008; Ramsayer et al. 2011; Kaltz et al 2012), and 
numerous models have been proposed to explain TL 
under various assumptions (e.g. Anderson et al. 1982; 
Ballantyne IV 2005; Engen et al. 2008). However, there is 
no consensus about why TL is so widely observed, how 
its estimated parameters should be interpreted in terms 
of underlying population dynamics, and when it might 
fail to be valid. 

Density-mass allometry is seen in at least two different 
contexts. First, the allometry is widely observed across 
taxa, when each taxon is described by an average body 
size and a population density, although the exponent of 
the power law appears to differ for different groups of 
organisms (Damuth 1981, 1987; Lawton 1989; Marquet 
et al. 1990; Silva and Downing 1995; Dunham and Vin- 
yard 1997; Enquist et al. 1998; Hendriks 1999; Schmid 
et al. 2000; Morand and Poulin 2002; Niklas et al. 2003; 
Makarieva et al. 2005; Reuman et al. 2008, 2009). This 
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form of density-mass allometry is sometimes referred to 
as Damuth's law. Second, density-mass allometry is often 
observed using the densities of individuals grouped by 
body mass irrespective of taxon. In plant ecology, this 
form of density-mass allometry is called the self-thinning 
law (Adler 1996); in marine ecosystems, it is called a size 
spectrum (Sheldon and Parsons 1967), and this is the 
form of density-mass allometry used in this article. 

Although the allometry of density with taxon average 
body mass (Damuth's law) has been controversial 
(Marquet et al. 1995, 2005), the allometry of density with 
individual body mass in marine ecosystems is clearly sup- 
ported by empirical data showing that the total biomass in 
logarithmic bins of body mass is approximately the same 
across a wide range of body sizes (Sheldon et al. 1972, 1977; 
Boudreau and Dickie 1992; Kerr and Dickie 2001). Equiva- 
lently, the density of organisms of a given mass (per unit 
volume per unit body mass) is a power-law function of 
body mass with an exponent close to —2 (Sheldon and Par- 
sons 1967; Piatt and Denman 1978; San Martin et al. 2006). 

Dynamic models of size spectra in marine ecosystems 
are based on a size-specific account of predation and 
growth as organisms eat one another, typically using the 
McKendrick-von Foerster partial differential equation 
(Benoit and Rochet 2004; Andersen and Beyer 2006; Blan- 
chard et al. 2009; Law et al. 2009; Hartvig et al. 2011; 
Zhang et al. 2012). Such models are deterministic and 
can display at least two modes of behavior, depending on 
parameter values. In the stable mode, the biomass in any 
bin of body size converges to a fixed limit in time. In the 
oscillatory mode, the biomass in any bin of body size 
converges to a periodic cycle. Other dynamics, such as 
divergence or chaos, may be possible and have not yet 
been investigated in detail. 

Because the oscillatory mode of dynamic size spectrum 
models predicts temporally varying population densities, 
it is natural to ask, for each body-size bin, (a) how the 
variance of population density is related to the mean of 
population density (does the model obey TL?), (b) how 
the mean population density is related to body size (does 
the model obey density-mass allometry?), and (c) how 
the variance of population density is related to body size 
(does the model obey variance-mass allometry?). Question 
(b) has been investigated in depth (Sheldon et al. 1972, 
1977; Boudreau and Dickie 1992; Andersen and Beyer 
2006), but neither question (a) nor question (c) has. 
Positive answers to these questions would provide novel 
interpretations of TL and variance-mass allometry in 
terms of the dynamic processes in size spectrum models. 
Negative answers would challenge the application of TL 
to dynamic size spectrum models. 

The relationships among variability, density, and body 
size in marine ecosystems are important from a practi- 



cal perspective because fish stocks appear to show 
increased variation under pressure of fishing (Hsieh 
et al. 2006; Andersen and Pedersen 2010; Rochet and 
Benoit 2012). However, the natural scaling of variability 
with density and body size in these ecosystems is largely 
unknown. Moreover, the traditional method of manag- 
ing fisheries, protecting young, small fish and harvesting 
old, large ones is coming under increasing scrutiny. One 
suggestion is that moving away from highly selective 
fishing, and instead spreading the fishing effort widely 
over species and sizes, would ensure a sustainable fish- 
ery while reducing waste and conserving biodiversity. 
This approach to fishing is referred to as balanced 
harvesting (Zhou et al. 2010; Garcia et al. 2012). Theo- 
retical results from a size spectrum model suggest that 
this approach has the potential to reduce the disruption 
to the natural size structure of the system (Law et al. 
2012), but the effects of harvest patterns on the scaling 
of variability with body size and density remain 
unknown. 

Here, we examined the extent to which TL and vari- 
ance-mass allometry apply to size spectrum models. We 
did this first using a mathematical argument about 
dynamics close to the boundary between stability and 
instability. Second, we used numerical methods applied to 
a model given by Law et al. (2012) over wide ranges of 
assumptions about predator-prey mass ratios, life histo- 
ries, and the forms and intensity of fishing. This model is 
based on an extension of the McKendrick-von Foerster 
equation to include a diffusion term (Datta et al. 2010) 
and a model for size-dependent reproduction (Hartvig 
et al. 2011). Our results showed that the model of Law 
et al. (2012) robustly predicted that size spectra should 
conform to TL, density-mass allometry, and variance- 
mass allometry. Moreover, simulation of a size-at-entry 
fishery, in which only fish above a minimum body mass 
can legally be caught, disrupted these relationships more 
than did balanced harvesting. 

Methods 

Size spectrum models are known to have equilibrium 
(steady state) solutions Us(x) and, in some cases, periodic 
solutions u{x, t) (Benoit and Rochet 2004; Andersen and 
Beyer 2006; Datta et al 2011). The variable u represents 
the density (per unit volume of water) of organisms of 
body mass w = Woe^, where Wq is an arbitrary mass, say 
that of an egg. Throughout u is thought of as a function 
of the logarithm x of body mass w. 

The solution of interest in this article was the periodic 
one, as this provided variation in density over time. 
We defined the mean M(x) and variance V{x) of the non- 
dimensionalized density at body size x over time as: 
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M{x) 



UopT 



u{x, t)dt 



V{x) = I u{x, tydt - M[xY 



Here, Mq^ was a reference density corresponding to the 
density of plankton at fish-egg size (Table 1). In the equi- 
librium case, M(x) equaled the equilibrium density U;{x) 
and V(x) was zero. 

Where the model had periodic behavior, but was close 
to the bifurcation from a stable equilibrium, analytical 
methods were used to predict the exponent for TL in an 
idealized form of the model (Datta et al. 2011; Plank and 
Law 2011). To investigate what happened in the more 
realistic model of Law et al. (2012), and when the system 
was not close to the bifurcation point, numerical methods 
were used to solve the model until it converged either to 
a stable equilibrium or to a periodic solution. In the peri- 
odic case, convergence was tested by numerically integrat- 



Table 1. Model parameters and values. 



Model parameter 



Symbol Value 



Dynamic size spectrum 

Biomass conversion efficiency 
Searcin rate scaling exponent 
Search rate constant 
Preferred log predator-to-prey 
mass ratio 
Diet breadth 

Intrinsic mortality constant 
Intrinsic mortality scaling 
exponent 

Mass at onset of senescent 
mortality 

Senescent mortality scaling 
exponent 
Fish life history 
Mass of egg 

Mass at maturation midpoint 

Asymptotic mass 

Reproduction function exponent 

Controls the width of the 

transition from 

immaturity to maturity 
Fixed plankton size spectrum 

Maximum body mass of plankton 

Plankton density at Xq 

Plankton size spectrum exponent 
Fishing 

Fishing mortality 

Minimum capture size 



6 



0.2 
0.8 

600 m^ year^'g" 
5 

1.5 

1 year"' 
0.25 



Wo exp(Xs) 440 g 



Wo 

Wo exp(x„ 



0.001 g 
150 g 



Wo exp(x„) 1 000 g 



P 

Pm 



0.2 
10 



Wo exp(Xp) 0.001 g 
Uqp 100 m"^ 



F 0-2.5 year" 

Wo exp(Xf) 1 50 g 



ing the system for a sufficient period of time and 
checking for periodic behavior of the total biomass of the 
system. The mean M(x) and variance V{x) of the periodic 
solution were then calculated using the above expressions. 
These expressions would be exact if the length of the inte- 
gration r were an exact multiple of the period of the 
solution. However, the period was not known a priori 
and was difficult to determine accurately. Hence, in 
numerical calculations, an approximate estimate for the 
period was obtained and then T was set to be 10 times 
this estimate. The results were also checked for conver- 
gence with respect to T. 

A linear relationship was fitted between the logarithm 
of the mean density In M(x) and the logarithm of body 
size x using standard linear regression: 

InM(x) = Oo + aix density-mass allometry (1) 

Here, x = 0 corresponds to Wq = 0.001 g (Table 1). 
Where the attractor was a periodic solution, relationships 
between the variance of density and body size, and 
between the variance and mean density were also fitted 
by the same method: 

In V{x) = bo + bix variance-mass allometry (2) 



In Vix) = Co + Ci\nM(x) Taylor's law 



(3) 



The regression parameters a^, bg, and Cq are referred to 
as intercepts. The parameters Oi, b^, and Ci correspond to 
the exponents of density-mass allometry, variance-mass 
allometry, and TL, respectively. The body size variable x 
spanned the range 0 < x < (Xj = 13) in steps of size 
Sx = 0.1 giving « = 130 data points for each linear regres- 
sion. The same range for x was used for each combina- 
tion of parameter values investigated. 

For the parameter values used by Law et al. (2012), the 
size spectrum converged to a stable equilibrium. Increasing 
the mean predator-to-prey mass ratio (PPMR) /? or decreas- 
ing the diet breadth a generally moved the model into 
regions of parameter space for which the equilibrium was 
unstable and the attractor was periodic (Plank and Law 
2011). We adopted a smaller value for the diet breadth a 
than that used by Law et al. (2012) to focus attention on 
periodic solutions. AU other parameter values were the same 
as those used by Law et al. (2012) unless otherwise stated 
(see Table 1). The life history parameters in Table 1 corre- 
spond to a species that eats both plankton and smaller fish, 
has an egg mass of wq = 1 mg, a maturity ogive with a mid- 
point at 150 g, and an asymptotic mass of 1000 g. The size 
spectrum parameters in Table 1 have the following mean- 
ings. Twenty percent of prey biomass is converted into 
predator biomass (K = 0.2). Predator search rate is propor- 
tional to w"', where the exponent a = 0.8 is widely used in 
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size spectrum models (Benoit and Rochet 2004; Andersen 
and Beyer 2006; Blanchard et al. 2009) and is based on cal- 
culations of how cruising speed of fish scales with body 
mass (Ware 1978). Predators typically consume prey items 
whose body mass is between 1/10 and 1/1000 of their own 
mass (/? = 5, a = 1.5). The largest source of mortality is 
predation {A = 600 m^ year^^ g ")> but there is in addi- 
tion an intrinsic mortality rate that is proportional to w^"'^^ 
(Brown et al. 2004) and a senescent mortality rate that 
applies to fish larger than 440 g. The senescent mortality 
rate is needed in a single-species model to prevent the 
buildup of large organisms that have no predators, but may 
not be needed in a community model containing other spe- 
cies that grow to larger body masses and continue to act as 
predators (e.g. Andersen and Beyer 2006). See Law et al. 
(2012) Appendix for full details of the model equations and 
parameters. 

Results 

In an idealized form of the model, designed to be mathe- 
matically tractable (Datta et al. 2011), the amplitude of 
perturbations to the equilibrium was a fixed proportion 

(a) = -4.44, a, = -0.71 , = 0.89 

Ol ■ ■ ■ 1 




-15 



0 5 10 15 

X 

(c) c„ = -2.3, c, = 1 .78, = 0.99 
-10| • • • 



of the mean density, near the bifurcation between equilib- 
rium and periodic behavior. If the ratio of perturbation 
amplitude to mean density (which is proportional to the 
coefficient of variation) were exactly constant, then the 
variance would be exactly proportional to the mean 
squared, which is TL with an exponent of exactly 2. The 
results of Plank and Law (2011) showed that the ratio of 
perturbation amplitude to mean density was not constant, 
but increased slowly with body size, and therefore 
decreased slowly with mean density. This predicted that 
the size spectrum should obey TL with an exponent of 
1.87, slightly <2 (see Appendix A for details). To obtain 
this analytical result, it was necessary to assume an infi- 
nite range of body masses, and reproduction was not 
explicitly included. However, the result provided a formal 
basis for the numerical calculations that follow. 

In numerical solutions of the model of Law et al. 
(2012), the three fitted linear regressions (1), (2), and (3) 
for the parameter values given in Table 1 all had coeffi- 
cient of determination > 0.88, indicating that the 
model solutions fitted a linear regression reasonably well 
(Fig. 1). The value of the TL exponent Cj was 1.78, which 
is close to the theoretically predicted value of 1.87. If all 

(b) fe„ = -10.2, fc, =-1.26, = 0.88 



-10 




0 5 10 15 

X 



-15 




-25 



-15 -10 -5 0 

\n{M) 

Figure 1. The model solution (solid curve) and the fitted linear regression (dashed straight line) for: (a) mean population density against log body 
size (x) (or density-mass allometry (1)); (b) variance of population density against log body size (or variance-mass allometry (2)); (c) variance of 
population density against mean population density (or Taylor's law (3)). Parameter values were as in Table 1. 
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three power laws held exactly, the parameters of density- 
mass aUometry, variance-mass allometry, and TL would 
be related exactly by = UiCi and = Cq + CiUq. The fit- 
ted values of b^ and foi agreed with the values of b^ and 
bi estimated using these relationships to the number of 
significant figures given in Fig. 1. 

To check how sensitive these results were to model 
parameters, we varied several key parameters within 
ranges that gave periodic behavior. For each set of param- 
eter values, we calculated the regression parameters and 
the value of r' for (1), (2), and (3). For example, when 
we varied the diet breadth a, the fitted exponents were 
fairly insensitive to the value of a (Fig. 2a). The intercept 
of density-mass allometry (1) was also fairly insensitive. 
However, the intercepts of variance-mass allometry (2) 
and TL (3) were more sensitive (Fig. 2b) because the 
model underwent a bifurcation at cr w 1.59. When the 
parameter a went above the bifurcation point, the system 
converged to a stable equilibrium rather than a periodic 
solution. The variance was therefore zero in this region of 
parameter space, and the curves (2) and (3) representing 
variance of population density were no longer defined, 
hence absent from the graph. As a approached the bifur- 



cation point from below, the intercepts of (2) and (3) 
tended to — oo. The values of (Fig. 2c) showed that the 
model solutions fitted TL (3) very closely (r^ > 0.98) for 
all values of a tested. Relationships (1) and (2) fitted less 
well, particularly for lower values of a but still had 

> 0.74 for all values tested. The fitted regression 
parameters agreed well with the parameters estimated via 
bi = fliCi and bo = Ca + CiUq across the range of values of 
a shown in Fig. 2. 

We repeated this analysis using several different model 
parameters: the preferred log PPMR (/?); the maturation 
midpoint (x^a,); the maximum size of the fixed plankton 
spectrum (xp); and the slope of the fixed plankton spectrum 
()')■ The maturation midpoint is the body mass at 
which 50% of fish have reached reproductive maturity. 
When we varied this parameter, we also varied the asymp- 
totic (maximum) body mass (x„o) so that the ratio of 
asymptotic mass to mass at maturation was constant. Some 
theoretical and empirical evidence suggests that this ratio 
varies little across species (Beverton 1992; Charnov 1993). 
Varying each of these parameters changed the linear regres- 
sion statistics (Table 2) remarkably little and the linear 
regressions fitted well > 0.7 for the relationships with 




0.6 



0.4 
0.2 



1.3 1.4 1.5 1.6 

Diet breadth (a) 

Figure 2. Effect of varying tine diet breadth a on the fitted linear regressions (circles - density-mass allometry (1); crosses - variance-mass 
allometry (2); squares - Taylor's law (3)): (a) exponents Si, fa,, and Ci; (b) intercepts ao, bo, and Cq; (c) r^. All other parameter values were as in 
Table 1. Dashed vertical line shows the value of a = 1.5 corresponding to the results shown in Figure 1. 



3172 



2012 The Authors. Published by Blackwell Publishing Ltd. 



J. E. Cohen et al. 



Taylor's Law and Body Size 



Table 2. The effect of varying modei parameters on the fitted rela- 
tionships among body size, mean population density, and variance in 
population density. 



Parameter Mean: 


Variance: 


Variance: 


and range body size (1) 


body size (2) 


mean (3) 


Diet breadth a [1.3, 1.65] 






Exponent [-0.7, -0.6] 


[-1.4, -1.2] 


[1.7, 2] 


Intercept [-7, -4] 


[-00,-10] 


[-00, 2] 


1^ [0.73, 0.89] 


[0.76, 0.91] 


[0.98, 1.00] 


Preferred predator-to-prey mass ratio ji [4, 6] 




Exponent [-0.8, -0.5] 


[-1.4, -1.1] 


[1.6, 2.1] 


Intercept [-7, -4] 


[-»,-9] 


[-00, 3] 


1^ [0.79, 0.90] 


[0.81, 0.89] 


[0.97, 1.00] 


Maturation midpoint [10.5, 13] 






Exponent [-0.8, -0.5] 


[-1.4, -1.2] 


[1.7, 1.9] 


Intercept [-5, -4] 


[-00,-10] 


[-00,-2] 


[0.85, 0.98] 


[0.85, 0.92] 


[0.98, 1.00] 


Maximum plankton size Xp [0, 3] 






Exponent [-0.8, -0.6] 


[-1.3, -1.2] 


[1.6, 1.8] 


Intercept [-5, -3] 


[-», -9] 


[-00,-2] 


1^ [0.88, 0.98] 


[0.88, 0.97] 


[0.96, 1.00] 


Plankton spectrum exponent [- 


2.2, -1.5] 




Exponent [-0.8, -0.6] 


[-1.4, -1.2] 


[1.7, 1.8] 


Intercept [-6, -4] 


[-14, -9] 


[-4, -2] 


[0.86, 0.90] 


[0.86, 0.88] 


[0.99, 1.00] 



Lower bounds of -oo for intercepts indicate that the variance tended 
to zero within the investigated parameter range. 



body mass; r' > 0.96 for TL) in all cases. For many of the 
parameters investigated, the system underwent a bifurca- 
tion from periodic to equilibrium behavior at some point 
in the parameter range, giving intercepts that tended to — oo 
at the bifurcation point as observed above. 

Two different size-based fishing protocols affected 
these relationships differently. In a standard size-at-entry 
fishery, all fish below a certain minimum catch size are 
protected from fishing, in principle. We incorporated this 
into the model by assuming a constant fishing mortality 
rate F for all fish above this body mass. In a model of 
balanced harvesting, fish of a given body mass are har- 
vested in proportion to their natural productivity (Garcia 
et al. 2012). We adopted the same measure of productiv- 
ity as Law et al. (2012): productivity P{x) was defined to 
be the product of the somatic growth rate and the density 
of individuals of a given size at the unexploited equilib- 
rium. The fishing mortality rate at logarithmic body size 
X was then set as ^iy(x) = FP(x)/Pq, where f is a constant 
and Po is the productivity at the smallest harvested body 
mass, which was set to be 1 g. For each fishing scenario 
considered, for a range of levels of fishing intensity F, we 
also calculated the total time-averaged yield. This was the 
sustainable yield in the sense that it was calculated from 
the long-term state of the system. 



Under size-at-entry fishing, the system behaved period- 
ically across the range of fishing intensity F we investi- 
gated (Fig. 3). Under balanced harvesting (Fig. 4), the 
system bifurcated to stable equilibrium behavior as the 
fishing effort increased (variance-mass allometry intercept 
bo and TL intercept Cq tended to — oo in Fig. 4b). Neither 
fishing method had a large effect on the exponents of the 
allometric relationships investigated. Balanced harvesting 
did not have a substantial effect on (Fig. 4c), but size- 
at-entry fishing substantially reduced for density-mass 
allometry and variance-mass allometry. This fishing pro- 
tocol substantially altered the size spectrum away from 
the density-mass allometric power-law (1). The size-at- 
entry fishery had a greater effect on the intercept Aq of 
equation (1) (compare circles in Fig. 3b and Fig. 4b). A 
lower intercept and a comparable exponent indicated that 
size-at-entry fishing reduced the biomass of the remaining 
stock to a greater extent than did balanced harvesting. 
Finally, the yield for size-at-entry fishing (Fig. 3d) with 
these parameter values did not exceed around 
0.065 g m^^ year^\ In contrast, balanced harvesting gave 
yields of up to 0.85 g m^^ year^' (Fig. 4d). We have not 
optimized either fishing protocol to find the level of fish- 
ing effort or (in the case of size-at-entry) the minimum 
catch size that would give the maximum sustainable yield. 
Nevertheless, these results supported the findings of Law 
et al. (2012), which showed that balanced harvesting gave 
a higher maximum sustainable yield than size-at-entry 
fishing, and disrupted the ecosystem's natural size struc- 
ture much less. 

Discussion 

Equilibrium properties of size spectrum models have been 
the subject of much recent research (Andersen and Beyer 
2006; Andersen et al. 2008; Blanchard et al. 2009; Datta 
et al. 2010, 2011; Plank and Law 2011). The behavior of 
size spectra away from equilibrium has, in comparison, 
received less attention, although real systems often do not 
operate at equilibrium. Zhang et al. (2012) showed that 
incorporating species diversity (via the asymptotic body 
size trait) into the community model promoted stability 
and reduced the amplitude of oscillating solutions. Rochet 
and Benoit (2012) showed that size-selective fishing 
increases the amplitude of oscillations in the size spec- 
trum. 

Here, in one of the first studies to quantify the behav- 
ior of a size spectrum model away from equilibrium, we 
investigated allometric (power-law) relationships among 
body mass, mean population density, and variance of 
population density in a model of marine size spectra. We 
used mathematical arguments to predict that populations 
should conform to these allometric relationships in an 
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0 0.5 1 1.5 2 2.5 0 0.5 1 1.5 2 2.5 

F(1/yr) F(1/yr) 



Figure 3. Effects of size-at-entry fishing on the fitted linear regressions (circles - density-mass allometry, (1); crosses - variance-mass allometry, 
(2); squares - Taylor's law, (3)): (a) exponents a^, fai, and Ci; (b) intercepts Bq, bo, and Cq; (c) r^; (d) yield. Fishing mortality was zero below 
minimum catch size and was constant F above minimum catch size. In these graphs, minimum catch size was Xf = 11. 9, which is the same as the 
mean maturation mass. Other parameter values were as in Table 1 . 



idealized form of a size spectrum model. Within ecologi- 
cally realistic ranges for model parameters, numerical cal- 
culations confirmed that these relationships were 
remarkably robust, especially the relationship between 
variance and mean of population density (TL). Even 
when the allometric relationships of the mean and vari- 
ance of population density to body size were slightly 
weaker, TL stiU appeared to hold. The coefficient of 
determination (r^) for TL was >0.96 in all cases studied 
and the exponent was always between 1.6 and 2.1 and 
usually slightly <2. 

Some size spectrum models do not explicitly include 
reproduction, but assume a constant density of individu- 
als with body mass corresponding to the mass of an egg 
(Benoit and Rochet 2004; Blanchard et al. 2009; Law 
et al. 2009; Zhang et al. 2012). Other models explicitly 
couple the density of eggs to reproduction by mature 
individuals (Andersen et al. 2008; Hartvig et al. 2011; 
Rochet and Benoit 2012). Analytical results from a model 
without reproduction (Plank and Law 2011) predicted 
that the size spectrum should obey TL approximately. 



with an exponent slightly <2. However, the goodness-of- 
fit and the exponents of TL and variance-mass allometry 
could be assessed more accurately across a range of 
parameter values by avoiding the artificial constraint of 
zero variance in the density of eggs. In our numerical cal- 
culations, we therefore used a model that explicitly 
included reproduction (Law et al. 2012). Under this 
model, there is a natural interplay between variance in 
the density of eggs and the variance in the density of 
mature adults. 

Normally, TL applies to a set of populations, each with 
its own mean and variance of density. TL and variance- 
mass allometry have not before been demonstrated within 
a model of a single population. For instance, earlier work 
has shown that TL with an exponent 2 emerges in the dis- 
crete logistic growth model, Nf+i = rNf(l — Nf/K), when 
different populations have different carrying capacities 
K. If N(f) is a solution to the model with carrying capacity 
K then aN{t) is a solution to the model with carrying 
capacity aK. A straightforward consequence is that the 
variance must be proportional to the mean squared across 
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Figure 4. Effect of balanced harvesting on the fitted linear regressions (circles - density-mass allometry, (1); crosses - variance-mass allometry, 
(2); squares - Taylor's law, (3)); (a) exponents a,, b-,, and c,; (b) intercepts ao, bo, and Cq; (c) f^; (d) yield. Fishing mortality was proportional to 
natural productivity with constant of proportionality F. All parameter values were as in Table 1. 



a set of populations with different values of K (Ballantyne 
2005). Size spectrum models are different in that they dis- 
aggregate organisms by body mass allowing individuals to 
grow, often over several orders of magnitude of body 
mass, over the course of their lives as a result of eating 
smaller individuals. Hence, all the data on mean and vari- 
ance of abundance come from organisms of different body 
sizes within one population or community. The models 
we used did not assume a particular carrying capacity. The 
long-term biomass dynamics instead arose as a conse- 
quence of the interaction between the dynamic size spec- 
trum and its (assumed fixed) resource spectrum, as well as 
between different body sizes (predator and prey) within 
the dynamic spectrum. The finding that these dynamics 
conform to TL so closely across a wide range of parameter 
values was therefore new and surprising. 

Balanced harvesting, that is, spreading fishing effort 
across as wide a range of species and body sizes as possi- 
ble, has recently been proposed to be a more efficient and 
less harmful way of fishing than traditional fishing regula- 
tions such as size-at-entry (Zhou et al. 2010; Garcia et al. 
2012). In the context of a single-species size spectrum 
model, one model for balanced harvesting has been to 
match fishing effort with natural productivity across the 



range of body sizes in the population (Law et al. 2012). 
In practice, this means focusing more effort on the more 
productive, smaller fish and reducing fishing intensity on 
the larger fish. A key potential advantage of this balanced 
harvesting regime is that it is predicted to alter the rela- 
tionship between mean population density and body mass 
much less than more selective fishing methods (Law et al. 
2012). Our results added further weight to this conclusion 
by showing that balanced harvesting also disrupted the 
variance-mass allometry less than size-at-entry fishing. 
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Appendix A. Theoretical basis for 
Taylor's law in the size spectrum 
model 

The periodic behavior of the model arises via a Hopf 
bifurcation from stable equilibrium to stable periodic 
orbit. At the bifurcation, a complex conjugate pair of ei- 
genvalues /i has zero real part X = Mm. Close to the bifur- 
cation, the limit cycle solution can be expressed as a 
small perturbation to the steady state Us(x): 

u{x, t) = u,{x){l + Be{x)e'"" + Be{x)e-'"") (Al) 

where e{x) is the eigenfunction associated with the eigen- 
value A = ioj and the overbar denotes complex conjugate. 
5 is a complex constant whose magnitude grows in pro- 
portion toy/j/i — ^*|, where fi is the bifurcation para- 
meter and /(* is its bifurcation point (Strogatz 1994). 

From equation (Al), the mean population density is 
simply M{x) = u^x) and, using T to denote the period of 
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the limit cycle, the variance is 
1 /"^ 

V{x) =- u{x, tfdt - M{xf 
T Jo 

= - u,{xf / {BB{x)e'"" + Be{x)e-'""Ydt 
T Jo 

= -u,{xf / {Bh{xfe^""' + Bh{xfe-^''' 
T Jo 

+ 2BBe{x)e{x))dt 

= 2M{xf\B\^\F.{x)\^, 

Since 




df = 0 



for any non-zero integer value of n. Analytical expressions 
for the eigenfanction £{x) are available only in certain 
idealized forms of the model. In the special case where 
there is no reproduction, body sizes span an infinite 
range, and a constraint a = y — 1 is placed on the model 
parameters (a is the power-law exponent for the scaling 
of predator search rate with body size and y is the slope 
of the steady-state size spectrum), the eigenfunctions are 
plane waves of fixed amplitude (Datta et al. 2011): 

e{x) = e*^ = cos(fcx) + isin(foc), 

where /: is a real-valued constant. In this case, the magni- 
tude of the eigenfanction is constant |e(x)| = 1, which 
implies by equation (A2) that Taylor's law (TL) holds 
exactly with an exponent of 2. 

When the model parameters do not satisfy the 
constraint a = 7 — 1, numerical results (Plank and Law, 



2011) showed that, for parameter values compara- 
ble to those in Table 1, the magnitude of the eigenfunc- 
tion |e(x)| increases with x, approximately proportionally 
to e"^ where a = 0.07. The mean density M{x) is propor- 
tional to e*' ^ ''^ where y = 2.05, so 

|e(x)|ocM(x)"/<'-''> 

Hence, the results of Plank and Law (2011) combined 
with equation (A2) predict that the size spectrum should 
obey TL: 

V{x) = ClBfM(x)'+'"/<'-''' (A3) 

for some constant C. The predicted exponent for TL is 
2 + 2a/(l — y) = 1.87. The fitted exponent in the model 
of Law et al. (2012), which includes reproduction, 
(Fig. Ic) is 1.78. From equation (A3), the intercept Cq for 
TL (3) is 

Co = In C -|- In \B\^ 

Close to the bifurcation, the magnitude of B scales with 
— /J* (Strogatz 1994). Hence, the intercept for TL is 

Co = constant + In \ 

This explains the logarithmically shaped decay to —00 
of the intercept Cq as the parameter approaches the bifur- 
cation point in Fig. 2b and 4b. 

These arguments were based on an idealized form of 
the model that assumed an infinite range of body sizes 
and did not explicitly include reproduction. Nevertheless, 
they gave some analytical insight into why TL applied in 
a model of size spectra, and they provided qualitatively 
accurate predictions of the exponent and intercept. 
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